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ABSTRACT 


A  detailed  analysis  of  the  horizontal  boundary  layer  structure  of  the  warm  front  of 
an  open  ocean  explosive  cyclone  in  Intensive  Observation  Period  (lOP)  2  of  the  Exper¬ 
iment  on  Rapidly  Intensifying  Cyclones  in  the  Atlantic  (ERICA)  is  conducted.  Data  for 
this  study  consists  of  aircraft  data  averaged  over  one  minute  supplimentcd  by  satellite 
and  drifting  buoy  observations.  Analysis  of  surface  winds  and  fluxes  was  done  using  the 
Brown-Liu  Marine  PBL  model.  Results  show  a  PBL  which  differs  from  that  found  in 
typical  cyclones,  with  large  latent  heat  fluxes  south  of  the  warm  front  and  with  relatively 
weak  sensible  heat  fluxes  about  the  warm  front.  Boundary  layer  stratification  was  stable 
north  of  the  warm  front  and  unstable  south  of  the  warm  front  A  mechanism  for  moist 
frontogenesis  is  proposed  whereby  the  de.stabilizing  effects  of  the  latent  heat  flux  en¬ 
hances  frictional  convergence  along  the  warm  front.  These  fluxes  warm  and  moisten  the 
cyclone's  warm  sector,  enhancing  unstable  convection  along  the  warm  front  and  thereby 
enhancing  the  vertical  motion.  This  enhanced  vertical  motion  would  strengthen  the 
geostrophic  del'ormation  of  the  0,  gradient  and  potentially  enhance  cyclogenesis. 
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I.  INTRODUCTION. 


Maritime  explosive  cyclogenesis  poses  a  problem  for  mariner  and  meteorologist 
alike  with  its  unpredictable  sudden  severe  weather.  Explosive  cyclogenesis  is  defined  as 
a  central  pressure  drop  of  24  mb  in  24  hours  at  60°N,  geostrophically  adjusted  for  other 
latitudes  by  multiplying  the  deepening  rate  by  ^  (Sanders  and  Gyakum  1980). 

Forecasting  explosive  cyclogenesis  is  difficult  due  to  incomplete  understanding  of 
the  triggering  mechanisms  that  distinguish  it  from  non-explosive  development.  It  has 
been  shown  that  both  upper  level  and  boundaiy  layer  dynamics  contribute  to  explosive 
development  but  their  relative  importance  and  the  extent  to  which  they  interact  is  un¬ 
clear.  It  is  most  probable  that  explosive  cyclogencsis  has  no  single  trigger  but  is  rather 
the  result  of  one  of  many  possible  combinations  of  forcing  mechanisms.  The  upper  level 
inlluences  on  explosive  development  are  best  known  and  include  enhanced  cyclonic 
vorticity  advection  due  to  subsidence  induced  vortex  tube  stretching  (Boyle  and  Bosart 
1986)  and  enhanced  frontal  vertical  circulation  due  to  jet  streak  interaction  with  surface 
barociinicity  (L'cccllini  et.al.  1984). 

The  role  of  boundary  layer  processes  in  explosive  development  is  less  understood, 
but  preliminary  investigations  have  focused  on  surface  fluxes  of  heat  and  moisture  as 
forcing  mechanisms.  Sanders  and  Gyakum  (1980)  showed  that  explosive  cyclones  were 
winter  events  with  a  strong  preference  for  development  over  the  warm  waters  of  the 
Kuroshio  current  cast  of.lapan  and  the  Gulf  Stream  current  east  of  the  United  States. 
This  relationship  suggests  that  heat  and  moisture  fluxes  resulting  from  air-sea  temper¬ 
ature  difiercnccs  are  a  factor  in  explosive  development.  Observational  evidence  for  the 
influence  of  surface  heat  fluxes  on  explosive  cyclogenesis  was  provided  by  Davis  and 
Entanual  (1988)  in  their  investigation  on  the  propensity  of  the  Limited-area  Fine-mesh 
Model  (LFM)  to  undcrpredict  rapid  maritime  cyclogencsis.  Their  hypothesis  was  that 
the  model  failed  to  account  for  warming  due  to  surface  heat  and  moisture  fluxes.  Twenty 
nine  cases  of  explosive  storms  oil  the  east  coast  of  the  United  States  were  examined  and 
in  all  cases  the  dilference  between  predicted  and  actual  development  could  be  associated 
with  the  occurrence  of  warming  due  to  surface  fluxes.  A  concurrent  climatological  study 
of  the  same  region  from  1963-1977  showed  a  strong  correlation  between  regions  of  large 
sea-surface  fluxes  and  regions  of  explosive  cycTogenesis. 
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While  emphasizing  atmospheric  warming  due  to  surface  fluxes,  Davis  and  Emanuel 
(198S)  did  not  explicitly  define  the  flux  contribution  to  cyclogenesis.  Their  contention 
was  that  surface  heating  acts  to  raise  1000-500  mb  thickness  towards  a  potential  thick¬ 
ness  defined  by  the  sea-surface  temperature.  To  achieve  this  thickness  increase  the  sur¬ 
face  fluxes  must  directly  feed  the  ascending  air  of  a  cyclone.  Other  studies  suggest 
mechanisms  other  than  the  fluxes  are  responsible  for  this  atmospheric  warming.  Kuo 
and  Reed  (1988),  for  example,  concluded  that  surface  fluxes  played  a  minimal  role  in  one 
case  of  explosive  cyclogenesis.  Although  they  found  that  condensational  heating  asso¬ 
ciated  with  areas  of  hea\w  precipitation  led  to  rapid  development  through  enhanced  la¬ 
tent  heat  release,  the  surface  heat  and  moisture  fluxes  contributed  verx'  little  to  this 
latent  heat  release.  The  dilTering  impacts  of  the  flaxes  may  be  due  to  details  in  the  dis¬ 
tribution  of  surface  fluxes  and  boundar}'  layer  processes  in  developing  cyclones. 

The  distribution  of  surface  fluxes  in  west  coast  non-explosive  mid-latitude  storms 
was  examined  by  Fleagle  and  Nuss  (1985)  as  part  of  the  Storm  Transfer  and  Response 
Experiment  (STREX).  Surface  fluxes  were  found  to  have  relatively  large  positive  values 
in  the  cold  sector  behind  the  cold  front  but  small  or  negative  values  in  a  200  km  region 
of  the  warm  sector  in  advance  of  the  cold  front.  This  distribution  of  surface  fluxes  would 
act  to  oppose  cyclogenesis  by  reducing  the  baroclinicity  of  the  cyclone  through  the 
u  arming  of  the  cold  air  behind  the  cold  front. 

This  negative  effect  of  the  fluxes  on  cyclogenesis  differs  substantially  from  Bosart 
and  Lin's  ( 1984)  examination  of  the  1979  President's  Day  storm,  which  suggested  a  more 
positive  role  of  surface  fluxes  in  explosive  cyclogcnesis.  In  this  particular  case,  a  cyclone 
developed  exp]osi\  cl}  after  moving  offshore  over  the  Gulf  Stream  in  the  vicinity  of  Cape 
Hatteras.  They  concluded  that  surface  sensible  heat  fluxes  ahead  of  the  surface  low  acted 
to  pre-condition  the  environment  by  enhancing  the  low-level  baroclinic  zone.  Numerical 
model  studies  by  L’ccellini  et  al  (1987)  confirm  the  significant  role  of  fluxes  in  fostering 
explosive  development  of  this  cyclone.  The  distribution  of  surface  heating  and 
moistening  in  this  case  increased  low-level  frontogenesis  as  well  as  directly  feeding  the 
regions  of  latent  heat  release  in  the  mid-troposphere. 

Nuss  (1989)  examined  the  relationship  between  flux  distribution  and  cyclogenesis 
through  a  series  of  numerical  simulations  using  identical  initial  atmospheric  conditions 
but  diflerent  surface  fluxes  and  sea-surface  temperature  distributions.  The  control  case 
for  this  experiment  used  a  zonal  SS'f  distribution  with  upward  heat  and  moisture  fluxes 
to  the  west  of  the  surface  low  and  downward  fluxes  to  the  east,  which  is  similar  to  flux 
distributions  of  west  coast  systems.  A  simulation  using  a  sinusoidal  SS'f  distribution 


was  then  run  which  produced  additional  upward  fluxes  to  the  northeast  of  the  surface 
low.  similar  to  the  Presidents'  Day  storm.  A  comparison  of  simulation  results  after  24 
hours  showed  that  the  second  case  deepened  more  rapidly  than  the  control  (1.65  mb 
h'’  vice  1.47  mb  Ir').  The  control  case  showed  large  upward  fluxes  in  the  cold  sector 
behind  the  cold  front  and  near  zero  fluxes  in  the  warm  sector  in  front  of  the  cold  front. 
In  contrast,  the  second  case  had  large  upward  fluxes  not  only  in  the  cold  sector  but  also 
to  the  northeast  of  the  surface  low  in  advance  of  the  warm  front.  Additionally,  a  small 
upward  flux  was  found  throughout  the  warm  sector.  The  increased  deepening  rate  for 
the  second  case  was  attributed  to  the  destabilizing  eflects  of  upward  fluxes  of  heat  and 
moisture  on  the  boundary'  layer  which  resulted  in  increased  downward  momentum  fluxes 
to  the  north  of  the  warm  front.  This  in  turn  enhanced  the  associated  frontal  ageostrophic 
secondary  circulation,  leading  to  increased  frontal  convergence  and  vertical  circulations. 
Cyclogenesis  was  then  enhanced  through  an  associated  increase  in  the  release  of  latent 
heat. 

These  studies  indicate  that  the  planetary  boundaiy  layer  (PBL)  and  frontal  structure 
associated  with  cyclones  has  the  following  attributes: 

•  .A  strong,  stable  capping  inversion  above  a  well  mixed  boundary  layer  in  the  cold 
region  behind  the  cold  front  and  in  advance  of  the  warm  front; 

•  \\'cak.  or  neutral  stratification  in  tiie  warm  sector  and  unstable  stratification  in  the 
cool  sectors  of  a  cyclone; 

•  .A  narrow  updraft  of  warm  air  above  the  surface  cold  and  warm  front  position, 
where  there  is  no  well  defined  boundary  layer; 

•  .A  frontal  ageostrophic  secondary  circulation  with  the  surface  ageostrophic  wind 
moving  from  the  cold  air  to  the  warm,  with  upward  motion  in  the  warm  air. 

From  the  last  two  above  items,  it  is  clear  that  surface  flu.xes  must  be  in  the  warm  sector 
for  them  to  directly  feed  the  ascending  air  and  warm  the  atmosphere.  This  is  contrary 
to  the  study  by  Fleagle  and  Nuss  (1985)  of  west  coast,  non-explosive  storms,  which 
placed  the  largest  upward  fluxes  in  the  cold  sector  with  small  or  negative  surface  fluxes 
in  the  warm  sector.  However,  in  Nuss's  (1989)  simulation  of  the  elTects  of  surface  heat 
and  moisture  fluxes  on  cyclogenesis,  rapid  development  was  associated  with  upward 
surface  moisture  fluxes  in  the  warm  sector  and  boundary  layer  stability  differences  across 
the  warm  front  that  apparently  enhanced  frontal  convergence.  It  is  not  known  whether 
these  processes  occur  in  actual  explosive  cyclones. 

Investigation  into  explosive  cyclogcne'sis  has  been  hindered  by  a  lack  of  observations 
on  the  small  scale.  Most  studies  and  their  conclusions  are  based  on  synoptic  scale  anal- 
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ysis  with  limited  observations.  This  thesis  will  attempt  to  resolve  the  role  of  surface  heat 
and  moisture  lluxcs  and  horizontally  varying  boundars'  layer  structure  by  exaniining  the 
details  of  the  boundary  layer  structure  of  an  explosive  cyclone  in  the  warm  front  region. 

More  specifically,  this  thesis  has  the  following  objectives; 

1.  Map  the  detailed  horizontal  boundary  layer  structure  of  an  open  ocean  explosive 
cyclone  using  one  minute  averaged  aircraft  data. 

2.  Describe  the  distribution  of  surface  heat  and  moisture  fluxes  in  relation  to  the 
warm  front. 

3.  Describe  the  stability  structure  of  the  boundary  layer  across  the  warm  front. 

4.  Examine  the  eflect  that  the  distribution  of  surface  heat  and  moisture  Huxes.  and 
PBL  stability  structure,  have  on  the  frontogenesis  and  ageostrophic  secondary  cir¬ 
culation  of  the  warm  front. 
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II.  DATA  AND  MODEL  DESCRIPTION 


The  Experiment  on  Rapidly  Intensifying  Cyclones  over  the  Atlantic  (ERICA)  was 
devised  to  "determine  physical  mechanisms  and  processes,  and  their  critical  spatial  and 
temporal  combinations,  which  can  account  for  the  winter  time  phenomenon  of 
explosively  developing  over-ocean  atmospheric  storms."  (Hadlock  and  Kreitzbcrg  19SS). 
Conducted  from  1  December  19SS  to  1  March  1989,  the  area  of  interest  was  centered 
oITthe  east  coast  of  the  United  States  in  a  climatologically  favored  region  for  explosive 
cylogcnesis  extending  from  40°N  to  47°\  between  59‘'\V  and  69°VV  (Hadlock  and 
Kreitzbcrg  19SS).  The  ERIC.A  field  study  encompassed  six  explosive  cyclones  with  in¬ 
tensive  observation  periods  (lOPs)  of  aproximately  36  hours  covering  the  period  of  ex¬ 
plosive  development.  This  thesis  examines  a  portion  of  the  second  intense  observation 
period  (IOP-21  from  03UU  U'l  C  to  0700  E  TC  14  December,  which  was  the  initial  stage 
of  rapid  cyclogenesis. 

Tor  lOP-2  a  variety  of  data  were  obtained.  The  horizontal  PBE  structure  was  de¬ 
termined  by  a  NOA.A  WP-3D  which  was  tasked  to  document  the  PBL  s  ructure  by  Hying 
a  sawtoothed  flight  profile  across  the  warm  front.  Vertical  structure  was  to  be  measured 
by  aircraft  deployed  dropwindsondcs,  of  whicii  only  one  provided  good  data.  The  air¬ 
craft  data  were  supplemented  by  three  drifting  buoys  from  a  field  of  fifty  drifting  buoys 
at  2u0  km  spacing  which  measured  mean  sea-level  pressure,  air  temperature,  and  sea- 
surface  temperature  (SST)  on  an  hourly  basis.  Satellite  SST  measurements  and  standard 
ship  observations  were  also  available  and  used  in  this  study. 

The  lOP-2  aircraft  data  which  were  as  used  to  construct  detailed  horizontal  ana¬ 
lyses.  consisted  of  76  mcteorlogical  and  flight  parameters  recorded  at  one-second  inter¬ 
vals.  The  first  step  in  processing  the  data  for  this  thesis  was  to  time  average  the  data  over 
6()-second  intervals  in  order  to  reduce  the  data  set  to  a  more  managable  size  and  remove 
very  small  scale  variablity  that  could  not  be  resolved  in  the  horizontal  analysis.  Care 
was  taken  not  to  average  over  any  data  gaps  or  around  turns  where  wind  data  may  be 
questionable.  Horizontal  homogenity  was  maintained  by  rejecting  any  measurements 
taken  when  the  aircraft  was  climbing  or  descending  at  greater  than  5  m  S"’  or  if  the  air¬ 
craft  significantly  deviated  from  the  desired  altitude  of  350  m. 

.‘\  Cressman  analysis  scheme  was  then  employed  to  produce  gridded  data  fields  from 
these  processed  aircraft  observations.  The  Cressman  analysis  was  done  over  a  50(t  km 


by  350  km  grid  with  a  10  km  grid  spacing.  This  grid  relative  to  the  aircraft  track  is  shown 
in  Fig  1.  Five  analysis  passes  were  made  on  the  data  using  decreasing  scan  radii  of  2<>.0. 
15.0,  10. 0.  "^.0.  and  3.0  km.  This  thesis  treats  the  four-hour  data  acquisition  effort  as  a 
snapshot'  of  the  storm  and  no  corrections  for  system  translation  were  made.  .A  cor¬ 
rection  for  system  translation  would  entail  an  assumption  of  quasi-steady  state  that  vas 
not  justified  given  that  the  cylone  deepened  12  mb  over  the  four  hours  when  data  were 
taken.  It  is  believed  that  the  errors  associated  with  any  such  correction  may  be  com¬ 
parable  to  analysis  errors  with  no  correction.  While  this  assumption  about  the  error, 
when  no  translation  correction  is  made,  may  not  be  true  in  general,  several  factors 
minimize  the  errors  in  this  case. 

•  The  warm  front  was  quasi-staionary  during  the  flight  and  the  cyclone  moved  east¬ 
ward  along  the  warm  front.  Consequently,  structure  to  the  east  of  the  the  low  is 
not  influenced  by  substantial  translation  of  the  primary  feature  and  there  i*-  liti'e 
blending  of  substantially  ditlV.cnt  data,  which  gives  relatively  small  errors.  1  o  the 
wc'-t  of  the  low.  the  analysis  consists  of  pre-cyclone  and  post-cyclone  data  and  the 
error  is  larger  due  to  this  blending; 

•  I  he  aircraft  track  moved  to  the  east  with  the  system  so  that  the  majority  of  the 
data  are  in  the  warm  front  region.  In  addition,  the  translation  of  the  cyclone  was 
relatively  short  compared  to  the  data  coverage  to  the  east  of  the  cyclone.  This  fa- 
\  ors  the  analysis  of  the  w  arm  front  which  is  the  primary  focus  of  this  study. 

Consequently  the  analysis  in  the  warnt  frontul  region  are  bcliesed  to  be  representative 
of  the  structure  e\cn  w  ith  some  time  aveaging  of  this  structure. 

‘fhe  data  fields  were  then  passed  to  the  Brown-Liu  marine  planetarc  boundary  layer 
model  lor  calculation  of  surl'acc  winds  and  wind  stress,  and  surface  heat  and  moisture 
fluxes.  I  his  FBI.  model  is  a  two-layer  model  with  an  Fkman  Taylor  layer  over  a  loga¬ 
rithmic  surface  layer,  fhe  model  calculates  surface  fluxes  and  winds  given  an  upper 
boundary  condition  which  can  be  the  top  of  the  I’BL  or  within  the  Fkman  layer.  1  he 
flight  lesel  winds  are  taken  as  the  winds  within  the  Fkman  layer  from  which  surface 
w  inds  are  calculated,  fhe  air  temperature  field  is  required  for  thermal  wind  effects  while 
the  air-sea  temperature  difference  from  air  and  sea  temperature  fields  is  needed  to  cal¬ 
culate  FBL  stratification  and  surface  fluxes  (Brown  and  Liu  1982). 

fhe  Brown-Liu  model  requires  surface  air  temperature  and  relative  humidity  to  cal¬ 
culate  the  initial  surface  fluxes  to  determine  the  stability  and  iteratively  solve  for  the 
surface  fluxes.  This  required  that  the  flight  level  data  be  extrapolated  down  to  the  sur¬ 
face.  To  make  this  extrapolation,  a  well  mixed  PBL  with  constant  equivalent  potential 
temperature  (f),)  and  mtxing  ratio  fq)  was  assumed,  'fhe  temperature  (T)  can  then  be 
determined  b\  the  relationship 
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Figure  1.  K)P-2  Aiicinlt  track  relative  to  analysis  grid. 
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w  iici  c  p  is  tiic  pressure  at  lliglu  aititiKle.  I,  is  tlie  latent  heat  of  vaporisation,  ('/>  is  the 
cpceilic  heat  at  constant  pressure,  q  is  the  nrixing  ratio,  and  0,  is  equivalent  potential 
tcnipcrature  nKasured  at  the  aircrtilt.  To  determine  relative  humidity  (R)  at  tlic  surlacc, 
the  saturation  Najior  pressure  (c.)  was  calculated  using  the  derived  surl'ace  temperature 
field  and  the  following  relationship  between  vapor  pressure  and  temperature; 


7.SxT-27VK) 

(',  =  6.1078x10  i-S5S'i 

Saturated  mixing  ratio  (q,)  was  then  calculated  from  the  saturation  \apor  pressure  and 
surface  air  pressure  given  by  the  following  equation: 


q,  =  0.622( 
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Finally,  the  relative  humidity  field  was  calculated  using  the  aircraft  mixing  ratio  (q)  and 
the  derived  saturation  mixing  ratio  (q,)  using  the  standard  formulation  for  relative  hu¬ 
midity; 


To  determine  surfact  luxes, the  model  requires  that  the  air-surface  temperature  and 
the  aii-surface  humidity  diflerences  be  calculated.  The  necessary  SST  field  was  created 
using  multiple  data  sources  in  a  Cressman  analysis  scheme.  A  first  guess  SST  field  was 
derived  from  a  NFSDIS  Local-Scale  Analysis  Region  6  gridded  SST  field  with  a  14  km 
resolution.  Because  the  Tiros-N  series  polar  orbiting  satellites  only  retrieve  SST's  under 
dear  conditions,  the  age  of  the  observations  used  in  the  analysis  are  variable.  The  par¬ 
ticular  analysis  used  was  performed  at  0107  L’TC  14  December  19SS.  but  the  youngest 
satellite  data  are  at  least  [  2  hours  old.  and  the  majority  of  the  observations  were  132 
hours  old.  The  14  km  analysis  was  interpolated  onto  a  70  x  55  10  km  grid  centered  on 
the  50  X  35  10  km  grid  used  for  the  other  data  fields.  To  eliminate  any  boundary  effects, 
the  outer  10  grid  points  were  discarded  after  the  analysis  was  completed. 

Three  other  data  sources  were  used  in  the  analysis;  Radiometric  temperatures  from 
the  aircraft,  and  sea-surface  temperatures  from  ships  and  buoys.  The  aircraft's  down¬ 
ward  looking  radiometer  observations  were  then  screened  to  eliminate  erroneous  data: 

1.  All  readings  from  altitudes  greater  than  500  m  were  discarded  (water  vapor  and  or 

undetected  clouds  in  the  column  below  the  aircraft  can  give  a  cold  bias); 

2.  One-minute  averages  were  discarded  if  any  one  second  had 

a.  a  relative  humidity  greater  than  88%,  or 

b.  a  liquid  water  content  greater  than  0.3  g  kg"',  or 

c.  a  relative  humidity  greater  than  82%  and  a  liquid  water  content  greater  than 
0.2  g  kg  ' 

3.  Data  were  not  discarded  on  the  basis  of  aircraft  turns  or  ascent  descent  rate. 

Because  the  satellite  observations  were  largely  132  hours  old  and  also  because  the 
ship  and  buoy  data  outside  the  aircraft  flight  path  were  quite  sparse,  ship  and  buoy  data 
dating  back  to  oooo  L  l  <''  9  December  were  examined.  Ship  data  were  included  only  if 
they  were  w  ithin  approximately  3°  of  the  satellite  field  or  were  strongly  supported  by 
independent  aircraft  or  buoy  data.  Older  data  and  especially  old  data  near  strong  gradi- 
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ents  were  more  critically  examined.  Most  all  of  the  buoy  data  were  included,  except  in 
cases  where  redundant  older  data  were  discarded  in  favor  of  collocated  younger  data. 

The  final  analysis  was  performed  using  a  Cressman  scheme  with  a  limit  of  8  degrees 
on  the  difference  from  the  first  guess  imposed  as  a  final  check  on  the  data  quality.  Scan 
radii  of  10.0,  5.0,  2.5,  1.0,  and  0.5  km  were  used  on  five  successive  passes.  Given  the 
combination  of  sparse  and  clustered  observations,  these  scans  most  closely  replicated  a 
hand  analysis. 


9 


111.  IOP-2  SYNOPTIC  SUMMARY 


As  discussed  in  Challant  (1989),  10I’-2  is  characterized  by  the  development  of  se¬ 
veral  cyclones  associated  with  the  advance  of  seseral  short  wave  troughs  in  rapid  suc¬ 
cession.  1  he  first  cyclone  developed  well  to  the  south  and  underwent  a  brief  period  of 
moderate  deepening  before  it  weakened  rapidly.  While  the  first  cyclone  was  beginning 
to  weaken,  a  much  more  significant  cyclone  had  begun  to  develop  to  the  southwest  just 
off  the  North  Carolina,' V’irginia  coast.  Ihis  second  cyclone  underwent  very  rapid 
cyclogcnesis  to  become  the  dominant  feature.  I  he  0600  U’l  C  13  December  1988  (from 
here  on,  date  and  time  will  be  given  as  date  time  UTC.  e  g.  13  0600  Ul  C)  mean  sea-level 
pressure  and  isotherm  analysis  from  Sanders  (1989)  show  the  low-level  structure  associ¬ 
ated  with  the  first  cyclone  (Fig.  2).  The  cold  advection  evident  in  Fig.  2  is  associated  with 
the  992  mb  low  to  the  southeast  of  Newfoundland  at  43'N,  47°W  and  the  1028  mb  high 
off  the  New  Fncland  coast  at  48°N,  65°W. 


Figure  2.  lOP-2  13/0600  UTC  surface  analysis:  Sea-level  pressure  (solid)  in 
tenths  of  mb  where  the  leading  9  or  10  has  been  dropped  and  surface 
temperature  (dashed)  in  ’C.  (Sanders  1989). 
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By  13, 1800  U  rC.  both  surface  cyclones  are  evident  as  two  distinct  low  centers  em¬ 
bedded  in  a  large  area  of  cyclonic  circulation  (Fig.  3)  just  off  the  Fast  Coast.  1  he 
northern  cyclone,  with  a  central  pressure  of  1004  mb  is  located  in  a  region  of  very  little 
thermal  gradient  and  verv'  weak  warm  thermal  advection  east  of  the  low  (Fig.  3j  but  is 
apparently  near  the  north  wall  of  the  Gulf  Stream  where  a  substantial  surface  temper¬ 
ature  gradient  exists.  Din'erencing  the  isotherm  analysis  at  13,06  U'l  C  and  13,  IS  U  l  C 
reveals  a  region  of  greater  than  5°  C  warming  over  12  hours  to  the  north  and  east  of  the 
cyclone.  1  he  ma.ximum  warming  occurs  just  east  of  the  northern  low  pressure  center 
at  about  37.5°N,  which  is  near  the  location  of  the  Gulf  Stream.  Although  this 

analysis  is  only  approximate,  this  rapid  warming  in  a  region  of  weak  thermal  advection 
and  near  the  the  maximum  SS  f  gradient  is  certainly  suggestive  of  warming  from  the 
surface  fluxes  associated  with  the  Gulf  Stream.  Substantial  upper-level  vorticity 
advection  is  found  over  this  northern  cyclone  in  synoptic  scale  500  mb  analyses  (not 
shown).  Hence,  there  is  ample  support  for  cyclogenesis  to  occur  in  this  region. 
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At  14,0000  L' TC,  this  northern  low  at  37'’N,  70°W  with  a  central  pressure  of  996 
mb.  began  explosive  development.  Examining  the  14,0000  U  I  C  mean  sea-level  pressure 
and  isotherm  analysis  (Fig.  4),  warm  advection  has  increased  at  the  surface  just  cast  of 
this  cyclone.  However,  the  largest  surface  warming  is  found  to  the  northeast  (40’N, 
64’W)  where  little  warm  advection  is  found.  Extremely  rapid  warming  of  5°  C  in  6  hours 
in  the  region  to  the  nortlieast  is  highly  suggestive  of  heating  due  to  surface  fluxes.  A 
more  substantial  thermal  gradient  develops  as  the  cyclone  intensifies  which  may  be  due 
in  part  to  surface  warming  by  the  heat  Ilux.  I  he  position  of  the  low  center  coincides 
with  the  northern  boundary  of  the  Gulf  Stream  as  indicated  by  satellite  SS T  analysis. 

1  his  relationship  between  the  surface  low  center  and  the  Gulf  Stream  was  maintained 
during  most  of  the  explosive  development  which  yielded  a  cyclone  track  along  the  Gulf 
Stream  (Eig.  5).  It  is  possible  that  such  movement  would  enhance  cyclogcnesis  as  the 
increased  surface  fluxes  over  the  Gulf  Stream  increases  the  frontogcntical  forcing  of  the 
agcostrophic  secondary  circulation  associated  with  the  warm  front.  Detailed  analysis 
using  the  aircraft  data  will  explore  this  question  in  the  next  section. 


Eiguie  4,  IOP-2  14/0U0()  UTC  surface  analysis:  Same  as  figure  2  (Sanders  1989). 
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During  the  next  12  hours  this  cyclone  deepened  30  mb  to  become  a  968  mb  low  at 
14,  1200  U  rC  as  sliown  in  Fig.  6.  I  he  first  cyclone  to  the  south  stopped  developing  and 
was  rapidly  absorbed  into  this  much  more  intense  northern  cyclone.  'Ihe  low-level 
structure  at  this  time  is  typical  of  many  mature  cyclones  with  rvell  defined  warm  and  cold 
IVonts  and  the  beginning  of  the  occlusion  process.  1  he  cyclone  fills  during  subsequent 
time  periods. 


Figure  5.  IOP-2  Cyclone  track:  SS  I  isotherms  in  'C. 


To  summarize,  lOP-2  was  a  multiple  low  event  triggered  by  the  olfshore  movement 
of  two  upper  air  short  wave  troughs  12  hours  apart.  A  surface  low  developed  cast  of 
I  lorida  in  response  to  the  first  trough.  'Ihe  second  upper  air  trough  triggered  the  for¬ 
mation  of  two  new  lows  in  a  surface  trough  that  extended  north  of  the  initial  system. 

1  he  northernmost  system  developed  explosively  w  hile  moving  along  the  northern 
boundary  of  the  Gulf  Stream  in  a  region  that  had  experienced  rapid  w^arming  over  the 
prccceding  18  hours,  wl)i!c  tlie  other  two  systems  weakened  and  were  assimilated  into 
this  major  cyclone.  Crude  warming  charts  constructed  using  the  temperature  change 
over  successive  isotherm  analyses  suggest  that  surface  (luxes  along  the  Ciuif  Stream  po¬ 
tentially  contributed  to  this  extraordinary  cyclone  and  will  he  explored  in  detail  in  the 
follow  rig  section. 
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rigiiie  6.  I()r-2  14/1201)  Ul  C  surface  analysis:  Same  as  figure  2  (Sanders  1989). 


IV.  BOUNDARY  LAYER  ANALYSIS 


From  the  synoptic  examination  in  the  previous  section,  the  pattern  which  emerges 
is  that  the  cold  advection  associated  with  the  preceeding  cyclone  produced  large  surface 
heat  and  moisture  Duxes  that  feed  the  developing  cyclone  in  lOP-2.  The  Duxes  may  have 
contributed  to  warm  frontogenesis  along  the  Gulf  Stream  where  the  cyclone  tracked. 
The  aim  of  this  section  is  to  examine  in  detail  the  surface  Dux  distribution  and  boundary 
layer  structure  around  the  warm  front  during  the  initial  phase  of  explosive  development. 
As  described  in  section  II.  the  aircraft  data  were  analyzed  on  a  grid  for  this  analysis. 

The  detailed  surface  structure  at  around  14  0600  UTC  is  shown  in  Fig.  7.  The  sur¬ 
face  low  and  frontal  positions  are  taken  from  a  separate  hand  analysis  of  the  aircraft 
observations  of  winds  and  pressure.  The  horizontal  surface  air  temperature  distributioji 
along  the  warm  front  shows  a  3°C  gradient  from  I8“to  21'’C  approximately  coincident 
with  the  SST  gradient.  Also  apparent  is  the  cold  front  extending  north  south  along 
66.5'’\V.  The  gridded  thermal  gradient  is  not  e.xactly  correspondent  to  the  frontal  anal¬ 
ysis  due  to  the  data  blending  resulting  from  system  translation  to  the  east  and  smoothing 
inherent  in  the  objective  analysis  method.  There  is,  however,  a  high  degree  of  confidence 
in  the  computed  distribution  about  the  warm  front  due  to  lack  of  blending  as  previously- 
detailed  in  Section  II.  V'ith  this  in  mind,  analysis  will  concentrate  on  the  warm  frontal 
and  warm  sector  structure  of  lOP-2. 

The  surface  heat  Dux  distribution  around  the  warm  front  (Fig.  8)  shows  downward 
Duxes  to  the  north  of  the  front  and  upward  Duxes  to  the  south.  The  magnitude  of  the 
downward  Duxes  to  the  north  of  the  warm  front  are  generally  small  (less  than  -50  W 
nr^)  with  more  significant  upward  Duxes  ( +  1()0  Wnr^)  to  the  south  of  the  front.  The 
small  magnitude  of  the  Duxes  north  of  the  front  indicates  that  the  air  over  the  cold  water 
is  in  near  equilbrium  with  the  sea-surface  temperature.  The  maximum  heating  essen¬ 
tially  coincides  with  the  SST  gradient,  which  suggests  the  advection  of  cold  air  over  the 
Gulf  Stream.  The  increase  in  Duxes  towards  the  east  suggests  that  the  air  is  being 
modilied  as  it  Dows  in  from  the  east.  This  Dux  distribution  is  opposite  to  the  distribution 
observed  by  Petersson  et.  al.  (1962)  and  Fleagle  and  Nuss  (1985)  where  cooling  occurred 
in  the  warm  air  and  warming  occurred  in  the  cold  air.  This  distribution  of  heating  results 
in  warm  frontogenesis  to  increase  the  baroclinicity  of  the  system,  which  is  similar  to  the 
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ill  "C  analyzed  from  aircraft  data.  Latitude  and  longitude  in  one  degree 
increments  with  a  center  latitude  and  longitude  of  38°N  and  66’\V. 


flux  influence  on  the  coastal  frontogcncsis  in  the  Presidents'  Day  storm  and  is  suggestive 
of  some  arctic  ice  edge  polar  lows. 

fhe  moisture  fluxes  (f  ig.  9)  are  considerably  stronger  than  the  heat  flux,  with  50U 
Wnr^  upward  fluxes  south  of  the  warm  front  and  relatively  weaker  upward  fluxes  to  the 
north,  fhe  magnitude  of  the  fluxes  in  the  cyclone's  warm  sector  are  very  large  and  ap¬ 
proach  values  recorded  in  cold  air  outbreak  studies.  However,  these  large  moisture 
fluxes  are  occurring  in  the  warm  sector  of  a  developing  cyclone  where  they  will  feed  the 
cloud  and  precipitation  region  of  the  cyclone.  As  with  the  heat  flux,  tlie  SS  f  and 
moisture  flux  gradients  are  essentially  coincident,  which  emphasizes  the  role  of  the  Gulf 
Stream  boundary  in  creating  the  flux  distribution.  In  addition,  the  moisture  flux  tends 
to  increase  towards  the  east,  which  suggests  that  the  easterly  winds  in  this  region  are 
bringing  in  relatively  dry  air  that  is  modified  by  the  moisture  flux. 

fhe  fact  that  these  strong  moisture  flu.xes  persisted  throughout  the  initial  phase  of 
explosive  development  implies  a  continuous  supply  of  cold,  dry  air  over  the  warm  watc'r. 
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Figure  8.  IOP-2  surface  sensible  heat  flux:  Sensible  heat  flux  in  W  m  ^  Grid  as 
in  Fig.  7. 

Such  a  situation  is  not  consistant  uith  the  expected  warm  ndvcction  from  the  south 
picdictcd  by  traditional  rnid-latitude  cyclone  theory.  Hxaniination  of  the  surface  wind 
anaivsis  (Fig,  |0)  reveals  strong  easterly  flow  with  a  wind  inaxiniuni  of  26  m  s  '  located 
south  of  the  warm  front.  1  hese  strong  easterly  winds  are  consistant  with  the 
fi ontogenesis  mechanism  proposed  by  Grotjalm  and  Wang  (1989)  whereby  air  from  the 
cold  sector  of  a  preceding  cyclone  is  drawn  into  the  warm  sector  of  a  following  cyclone. 
Ihe  l.FObOO  L  fC  synoptic  analysis  (Fig.  2)  shows  such  a  preceding  cyclone  to  the 
northeast  of  the  incipient  explosive  cyclone.  A  study  of  the  boundary  layer  over  a  larger 
region  that  includes  the  tracing  of  air  parcel  trajectories  would  be  required  to  confirm 
that  air  parcels  move  from  the  cold  sector  of  the  previous  low  into  this  region  of  the 
warm  front,  \\  hile  this  is  bes  ond  the  scope  of  this  study,  the  wind  flow  near  the  front 
in  the  small  scale  analyses  do  support  this  interpretation. 

As  discussed  above,  the  moisture  fluxes  arc  advantageously  placed  to  feed  the  warm 
frontal  updraft  region  and  thus  enhance  latent  heat  release.  As  suggested  by  Davis  and 
limanua!  ( 1989).  surface  heat  and  moisture  fluxes  may  act  to  destabilise  moist  symmet- 


17 


Figure  9.  IOP-2  surface  latent  lieat  lliix:  Latent  heat  (lux  (solid)  in  \V  m  ^  SS  I 
(dashed)  in  °C.  Grid  as  in  Fig  7. 

rically  neutral  regions  and  produce  deep  slantwise  convection  in  the  warm  Frontal  region. 
While  no  analysis  of  the  stability  was  done,  vertical  or  slantwise,  both  satellite  and  radar 
images  (not  shown)  show  strong  convection  along  the  warm  front  with  radar  returns  on 
the  order  of  50  db;7  being  recorded,  fhe  warming  and  moistening  by  the  fluxes  in  this 
region  may  contribute  to  the  convective  instability  and  clearly  supply  moisture  to  the 
clouds.  A  more  detailed  analysis  of  the  deeper  structure  is  required  to  quantify  the 
contribution  of  these  fluxes  to  the  convection. 

Nuss  (1989)  highlighted  the  f'rontogenetical  role  of  enhanced  horizontally  varying 
momentum  transfer  to  the  sea  surface  due  to  the  destabilizing  elfects  of  surface  fluxes 
in  the  boundary  layer.  While  conclusive  dropsonde  data  is  not  available  for  the  cyclone's 
warm  sector,  some  idea  of  boundary  layer  stability  can  be  inferred  from  comparing  flight 
level  winds  and  surface  winds.  South  of  the  warm  front,  surface  wind  speeds  are  equal 
to  flight  level  winds,  fhis  indicates  a  deep,  well  mixed  boundary  layer  with  flight  level 
winds  being  mixed  down  to  the  surface,  winch  agrees  with  the  upward  heat  fluxes  in  this 
region.  In  contrast,  surface  winds  north  of  the  warm  front  arc  only  one  half  the  strength 


18 


Figure  1(1.  IOP-2  suif.ice  Miiuls:  Isolaclis  in  5  m  s  '  intervals.  (Jrici  as  in  Fig  7. 


of  lliglu  level  winds,  indicating  a  more  stable  boundary  layer,  also  in  agreement  u  ith  the 
downward  Iluxcs  in  this  region.  I  he  less  stable  boundary  layer  in  the  cyclone  s  warm 
sector  would  enhance  the  transfer  o(  momentum  to  the  surlacc  as  suggested  by  Nuss 
(F)S9)  but  this  increased  momentum  trairsfer  is  on  the  opposite  side  ol’ the  warm  front 
as  found  in  the  numerical  simulations.  However,  the  cllcct  is  similar  in  that  the  resultant 
surface  winds  arc  more  convergent  than  those  at  the  Might  level  of  the  aircralt.  I  his 
frictional  convergence  is  Irontogcnetical  which  is  evident  from  tlie  observed  structure. 

1  he  boundary  layer  structure  that  emerges  from  this  analysis  near  the  warm  front 
dilfers  liom  that  described  in  west  coast  cyclones  by  I  Icagic  ct  al.  (19S8).  The  unstable, 
well-mixed  boundary  layer  occurs  in  tlic  warm  sector  and  is  apparently  maintained  by 
the  easterly  flow  of  dry,  cold  air.  A  more  stable  boundary  layer  occurs  north  of  tlic  warnr 
front  uhcrc  the  air  lias  apparently  reached  near  equilibrium  with  the  sea  surface.  Ibis 
structure  favors  frontogenesis  througli  direct  diabatic  licating  from  the  sea  surface  and 
through  its  influence  on  the  horizontally  varying  friction. 

I  0  further  highlight  the  potential  impact  of  this  boundary  layer  structure  on  the 
dynamics  of  the  cyclone,  a  qualitatise  examination  of  the  (^-vector  foim  of  the  omega 


equation,  as  given  by  Hoskins  et  al.  (1978)  was  made.  As  the  Q-vectors  will  be  used 
descriptively  and  no  vectors  are  actually  calcilated,  only  a  brief  discussion  of  Q-vectors 
will  be  provided.  A  complete  derivation  of  the  Q-vector  form  of  the  omega  equation  can 
be  found  in  Appendix  A.  The  Q-vcctor  form  of  the  omega  equation  is  written  as 


+ 


=  -2V.Q, 


where  o  =  is  the  static  stability  parameter  and 


Q  =  (e,  .(?;) 


This  form  of  the  omega  equation  ignores  the  effects  of  diabatic  processes  in  forcing 
vertical  motion.  11  the  atmosphere  is  saturated,  then  one  can  directly  include  the  elfccts 
of  latent  hear  release  on  tlie  \  crtica!  motion  by  considering  the  forcing  in  terms  of  moist 
baroclinic  dynamics.  Therefore,  a  moist  Q-vector  can  be  defined  by  assuming  that  in  a 
saturated  atmosphere  ^),  is  equisalent  to  0.  Substituting  Q,  for  6>  in  the  Q-vcctor  omega 
equation  results  in 

/  fi;  rl,  \ 

\  '  / 

whuh  gives  the  geostrophic  frontogentical  forcing  of  vertical  motion  in  a  saturated 
cm  ii  oment.  .A  more  complete  description  of  the  moist  frontogenesis  can  be  obtained  by 
examining  the  two  dimensional  frontogenesis  equation.  The  moist  frontogenesis 
equation  is 


V 


The  first  term  represents  non-conservative  moist  processes  and  is  zero  except  at  the  sea 
surface  where  surface  heat  and  moisture  fluxes  alter  0,  due  to  interaction  with  the 
boundary.  1  he  second  and  third  terms  represent  the  horizontal  geostrophic  deformation 
ol'  the  0,  gradient,  which  is  the  moist  Q-vcctor  forcing.  Although  the  geostrophic  forcing 
is  not  known  in  this  case,  the  total  wind  and  6,  distribution  suggest  moist  frontogentical 
forcing  by  the  horizontal  deformation  (Fig.  10).  In  the  warm  sector  south  of  the  warm 
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front,  large  (~500  W  iit^)  upward  surface  moisture  lluxes  and  upward  surface  heat  fluxes 
act  to  increase  surface  0,  and  strengthen  frontogcntical  forcing  of  the  warm  front,  fhis 
warming  and  moistening  of  the  warm  sector  also  destabilizes  the  boundary  layer  in  the 
cyclone's  warm  sector,  which  apparently  results  in  enhanced  convection  along  the  warm 
front.  This  enhanced  vertical  motion  would  then  strengthen  the  geostrophic  deformation 
as  indicated  by  the  Q-vector  omega  equation,  which  potentially  enhances  the 
cyclogenesis.  While  this  analysis  is  rather  speculative,  it  docs  suggest  the  possible  dy¬ 
namical  processes  by  which  the  surface  fluxes  enhance  cyclogenesis  in  this  case.  It  re- 
nuuns  to  be  shown  through  a  more  quantitative  analysis  whether  the  surface  fluxes 
contribute  as  significantly  as  suggested  by  this  qualitative  analysis. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS. 

The  detailed  boundary’  layer  analysis  of  lOP-2  has  shown  that  the  distribution  of 
surface  fluxes  about  the  warm  front  in  an  explosive  cyclone  provides  favorable 
frontogentical  forcing  to  increase  cyclogenesis.  Specifically,  this  study  has  shown: 

1.  The  horizontal  boundary  layer  structure  of  the  cyclone  during  its  explosive  phase, 
when  accurately  mapped  on  a  small  scale  using  aircraft  data,  has  a  boundary  layer 
structure  that  difi'ers  from  that  found  in  typical  cyclones. 

2.  Heat  fluxes  about  the  warm  front  were  relatively  small  due  to  air  and  sea  surface 
temperatures  being  in  near  equilibrium.  In  contrast,  large  upward  moisture  fluxes 
were  present  south  of  the  warm  front  with  weaker  fluxes  to  the  north  of  the  front. 

3.  Boundary  layer  stratification  was  stable  to  the  north  of  the  warm  front  and  unsta¬ 
ble  to  the  south  of  the  warm  front.  I  his  resulted  in  dilTerent  momentum  fluxes  on 
either  side  of  the  front  to  produce  convergence  in  the  surface  winds. 

4.  The  surface  heat  and  moisture  flux  distribution  as  well  as  the  wind  field  are  favor¬ 
able  for  moist  frontogenesis  along  the  SST  gradient.  Given  nearly  saturated  condi¬ 
tions.  this  enhanced  surface  frontogenesis  would  contribute  significantly  to  the 
cydogenesis. 

3Vhile  the  analysis  points  to  a  potentially  large  role  for  surface  fluxes  in  explosive 
cydogenesis.  the  lintitations  inherent  in  this  single  case  study  suggest  the  following  rec¬ 
ommendations  for  further  research; 

1.  Thi'i  anahsis  of  the  PBL  in  a  limited  region  should  be  extended  to  a  larger  domain 
and  include  a  trajectory  analysis. 

2.  .\  complete  .’-D  anahsis  should  be  done  to  examine  the  dynamics  that  are  driving 
the  agcostrophic  circulation. 

3.  The  time  evolution  of  the  cyclone  and  PBL  structure  should  be  determined. 

4.  Lxtend  this  analysis  to  LRICA  lOPs  with  intermittent  explosive  development.  .A 
comparison  of  the  surface  flux  distribution  and  magnitude  during  explosive  and 
non-explosive  phases  would  highlight  any  possible  differences  in  the  role  of  the 
boundar}  layer  processes. 
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APPENDIX  A.  Q-VECTOR  DERIVATION 

Derivation  of  the  Q-vector  form  of  the  omega  equation  begins  with  the  quasi- 
geostrophic  momentum,  thermodynamic  and  continuity  equations: 
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For  the  Q-vector  derivation  the  following  assumptions  are  then  made: 

1.  f  plane  =  0, 

2.  frictionless  F,  =  Fy  =  0, 

• 

3.  adiabatic  0  =  0  and 

4.  constant  static  stabilitv  -r—  =  constant. 

■  cp 

Momentum  equations  are  then  difTerentiatcd  bv  -7—  and  the  thermodvnamic 

r  r  ' 

equation  by  rry  or  For  the  x  component  equation: 

(^  +  =  0  (2) 
cy  \_\  cl  *  /  cp 

Equation  ( 1 )  becomes 
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/  3  -  \  du„  cVg  dv  g 

(4-+V  .v)-—  +  -^.vUg-f,-r^  =  0 

\  ct  ^  J  cp  cp  ^  cp 

and  equation  (2)  becomes 

{  c  ,7'  r7\  cd  ,  ,  do)  dd 

(  —  +  +  y— — •\d  +  y— - =  0 

\  Cl  °  J  cy  cy  cy  cp 

Expanding  the  second  term  of  equation  (3), 
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and  by  using  the  thermal  wind  relationships 


iH  =  ,sl.  .  zIl 

cp  ^  O'  ’  cp 


c9 


and  the  fact  that  the  horizontal  divergence  of  the  geostrophic  wind  is  zero,  (5)  reduces 
to: 

ci;  cr„ 

-  =  Q,. 

cp  -  cy 

1  his  term  will  be  defined  as  Qj  and  is  identical  to,  but  of  opposite  sign  to,  the  second 
term  of  (4).  (3)  and  (4)  can  then  be  written  as: 


(4  +  h.v),4^ -a  +  v-oi-P- =0.  (7) 

\  Cl  *  /  cy  cy  cp 

Since  the  first  term  of  (6)  and  (7)  are  equal  through  the  thermal  wind  equation,  sub¬ 
tracting  (7)  from  (6)  results  in 
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In  quasi-geostrophic  theory,  geostrophic  motion  tends  to  destroy  thermal  wind  balance 
which  ageostrophic  motion  acts  to  restore.  From  (8)  it  can  be  seen  that  Qi  represents 
a  restoring  force  on  geostrophic  balance  through  differential  ageostrophic  horizontal 
motion  acting  on  horizontal  winds  and  differential  vertical  motion  acting  on  the  hori¬ 
zontal  temperature  gradient. 

A  completely  analogous  derivation  for  the  y  component  produces 


-/o 


cp 


y 


c6  c(o 
cp  cx 


—2Q\  and 


(9) 


Qx 


oK 

-y  .  ve . 

cx 


To  obtain  the  omega  equation,  we  differentiate  (8)  by-^  and  (9)  by  rry  and  add  the 
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two  together: 
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.Adding  yields: 
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Continuity  is  then  used  to  obtain: 
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r'urthcr  simplification  results  in 
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where  o  = 


+  4-fr>  =  -Jv.Q, 
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-f^y  is  the  static  stability  parameter. 
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